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Abstract 


An NLO QCD analysis of the final HERMES data on pion multiplicities is 
presented and a new set of pion fragmentation functions is extracted from the best 
fit to the data. We have studied the so-called [x, z\ and z] presentations of their 
data, as given by HERMES, which, in principle, should simply be two different ways 
of presenting the experimental data. We have based our extraction on an excellent 
fit to the [Q^,z] presentation of the data. We also draw attention to what appears 
to be a problem with the [x, z\ presentation of the HERMES data. 


PACS numbers: 13.60.Hb, 12.38.-t, 14.20.Dh 


1 Introduction 


In the absence of charged cnrrent nentrino data, the experiments on polarized inclnsive 
deep inelastic lepton-nncleon scattering (DIS) yield information only on the snm of quark 
and antiquark parton densities (PDFs), Aq + Aq, and the polarized gluon density AG. 
In order to extract separately Aq and Aq other reactions are needed. One possibility is 
to use the polarized semi-inclusive lepton-nucleon processes (SIDIS) l + N^l' + h + X, 
where h is a detected hadron (pion, kaon, etc) in the hnal state. In these processes 
new physical quantities appear - the collinear fragmentation functions which 

describe the fragmentation of quarks and antiquarks into hadrons. Due to the different 
fragmentation of quarks and antiquarks, the polarized parton densities Aq and Aq can 
be determined separately from a combined QCD analysis of the data on inclusive and 
semi-inclusive asymmetries. The key role of the fragmentation functions for the correct 
determination of sea quark parton densities Aq, especially of the polarized strange quark 
densitiy, was discussed in [1] . Note that the W data from RHIC give no information about 
the polarized strange quark density and cannot help to solve the so called ’’strange quark 
polarization puzzle” (see the second reference in [1]). 

There are different sources to extract the fragmentation functions (FFs) themselves: 
semi-inclusive e~^ e~ annihilation data, single-inclusive production of a hadron h at a 
high transverse momentum in hadron-hadron collisions, unpolarized semi-inclusive DIS 
processes. It is important to mention that the data on hadron multiplicities in unpolarized 
SIDIS processes are crucial for a reliable determination of FFs, because only then can one 
separate Dq{z, Q^) from D^{z, Q^) (from the other processes only the sum of them can be 
determined). The hrst global analysis based on all these reactions was carried out by de 
Florian, Sassot, Stratmann (DSS) group [2]. As a result, the properties of the extracted 
set of FFs signihcantly differed, especially in the kaon sector, from those of the other then 
published sets of FFs |3] determined from analyses in which the SIDIS data have been not 
included. Unfortunately, the DSS FFs were based on the unpublished IIERMES’05 SIDIS 
data on hadron multiplicities |1], which were not conhrmed in the hnal HERMES data 
[S]. Indeed the hnal HERMES data diher signihcantly from those used in the analysis 
of [2] so that the FFs extracted in [2] are incorrect (see Figures 9 and 10 in [3] for LO 
DSS FFs, and hg. 5 in [9] for NLO DSS FFs). Moreover, in the extraction of the next- 
to-leading order (NLO) DSS set of FFs there was a mistake (see the correction in the 
Appendix in [6]) in the expression for the longitudinal gluon Wilson coefficient function 
in the theoretical formulae for the multiplicities (we, independently, became aware of this 
error recently). It has turned out that not only the DSS FFs, but all the other sets of pion 
and kaon FFs presented in [3] are NOT in agreement with the hnal HERMES [5] and the 
preliminary COMPASS data [7] on hadron multiplicities. In our paper [8] a theoretical 
analysis of these data was performed and new sets of pion fragmentation functions were 
extracted from the best NLO QCD hts to the data, and it was shown that they disagree 
signihcantly with the pion FFs determined from all previous analyses. Very recently de 
Florian at ah (DSEHS) have presented results on pion FFs obtained from their new global 
QCD analysis [9] using the hnal HERMES and the preliminary COMPASS data on pion 
multiplicities. 

In our paper [8] we pointed out a possible inconsistency between the HERMES [x, z] 
and presentations of their data on pion multiplicities. Bearing in mind that the 
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semi-inclusive DIS hadron production processes are essential for the separation of and 
D’i fragmentation functions, we present in this paper a more detailed discussion of our 
previous analysis in which we have taken into account the mistake in the longitudinal 
gluon Wilson coefficient function [10] present in our previous analysis, and have used the 
corrected version given in [Bj. Also, instead of the NLO MRST’02 set im we have here 
utilized the newer NLO MSTW’08 set [T2| of unpolarized parton densities and study the 
influence of this on the extracted FFs. 


2 QCD treatment of pion multiplicities 


The multiplicities z) of pions using a proton (deuteron) target are defined 

as the number of pions produced, normalized to the number of DIS events, and can be 
expressed in terms of the semi-inclusive cross-section and the inclusive cross-section 


Q^, z)/dxdQ‘^dz (x, Q^, z)/dxdQ'^dz 

d‘^N^^Ij^^{x,Q‘^)/dxdQ‘^ / dxdQ"^ 

(1 + (1 - yf)2xF^^^j^^{x, Q^, z) + 2(1 - y)xF^p(rf)(x, z) 

(1 + (1 - l/)^)2xFip(rf)(a:, + 2(1 - y)FLp(d}{x, 


In Eq. ([H) F[,Fl and Fi,Fl are the semi-inclusive and the usual nucleon structure 
functions, respectively. F[ and F^ are expressed in terms of the unpolarized parton 
densities and fragmentation functions (see 0 ). while Fi and Fl are given purely in terms 
of the unpolarized parton densities. 

We have assumed in our analysis that isospin SU{2) symmetry for the favored and 
unfavored fragmentation functions holds 


D:^{z,QI) = Df{z,Ql), Df{z,Ql) = Df{z,Ql), 


( 2 ) 


and in addition, the following relations for the fragmentation of strange quarks into a 
pion: 

Dfiz,Ql) = Dfiz,Ql) = Dfiz,Ql). (3) 

Due to the charge conjugation invariance of the strong interactions the fragmentation 
functions ^ can be expressed through 

= d;-(z,qI)^d;\z,q^„). (4) 

As a result, we have to extract only three independent FFs (Dj[^, from the 

NLO QCD fit to HERMES proton and deuteron data on pion multiplicities. The charm 
contribution to the multiplicities is not taken into account. In the theoretical analysis 
of the data the Mellin transform technique [10] was used to calculate the semi-inclusive 
Fl^j^{x,Q‘^,z) and the usual Fi^l^x^Q"^) nucleon structure functions in Eq.([TD from their 

moments. The expressions for the moments of the Wilson coefficient functions ^ (x, z) 
needed in these calculations can been found in HD]. As was mentioned in the Introduction, 
the error in the gluon Wilson coefficient, was corrected. Compared to our previous 
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fit | 8 ] where for the unpolarized PDFs we have used the NLO MRST’02 set [n] , we use 
now the NLO MSTW’08 set [12], for which the strange quark density s{x, is not equal 
to s{x,Q‘^). Note that we have chosen this set of PDFs in order to be able to compare 
correctly our extracted pion FFs with those of DSEHS obtained from the recent global 
ht [S] where the MSTW’08 set of PDFs has been used. The influence of the choice of the 
unpolarized densities on the extracted FFs will be discussed. 

For the input FFs the following parametrization sX Ql = 1 GeV^ was used: 


zDf{z,Ql) 


- zY^[l + 7j(l - z)Y 

B[ai + 1, + 1] + 7ji?[Q;j + 1, /3j + 5^ + 1] 


(5) 


where the parameters {W, /S*, 7 *, 6i} are free parameters to be determined from the ht 

to the data. Here, i stands for u,u and g, while B{a,b) denotes the Euler beta function, 
and the W are chosen in such a way that they represent the contribution of zD^ to the 
momentum sum rule. 


3 Results of analysis 

Let us discuss now our results on the pion FFs extracted from our NLO QCD ht to the 
HERMES proton and deuteron data on pion multiplicities, corrected for exclusive vector 
meson production [5]. In our study we have analyzed the [Q'^^z] and [x^z] presentations 
of the data (see Fig. 8 in [5], left column, second and third lines) for which the multiplic¬ 
ities do not depend on PhXi were Phi_ is the component of the hadron momentum, P/j, 
transverse to the momentum of the virtual photon. They correspond to two-dimensional 
projections obtained by the HERMES group from the full HERMES data sets [Q^, 2 ;, Phi\ 
and [x^z^PhYi respectively. The pion multiplicities are given for 4 z-bins [0.2-0.3; 0.3- 
0.4; 0.4-0. 6 ; 0.6-0. 8 ] as functions of the mean value of >, of each individual 

bin for the presentation or as functions of the mean value oi x, < x >, of 

each individual x bin for the [x^z] presentation. Note that for the [Q‘^,z] presentation 
there is no binning in x. This means that the multiplicity measured in a given bin, 
Qmin < < Qmax-: corresponds to the summing over all possible values of x belonging 

to the strip in the {x — QY plane, bounded by Qmax the kinematics of the 

HERMES experiment. And vice versa, for the [x, z] presentation, there is no binning in 
and the multiplicity measured for each x bin corresponds to all possible values of 
belonging to the {x, QY strip fixed by the boundaries of the x bin and the kinematics of 
the HERMES experiment. Thus, in principle, in the theoretical calculation of the pion 
multiplicities one has to integrate the semi-inclusive and inclusive cross section on RHS 
side of Eq. ([1]) over the x and regions corresponding to each bin for the [Q'^,z] 
presentation or to each x bin for the [x, z] one. It turns out however, that replacing x and 
by their mean values < x > and < > in the calculation of the multiplicities leads 

to very small difference. Further details are presented later. 

The total number of the 7r~^ and tt” data points for each of the presentations is 144, 72 
for 7 r+ and 72 for 7i~ data. In the case of [Q^, z] presentation of the data a good fit to the 
proton and deuteron data is achieved, y^/d.o.f = 123.95/132 = 0.94 for 144 experimental 
points and 12 free parameters. The errors used in the fit are quadratic combinations of 
the statistical and point-to-point systematic errors. We have found that the description 
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Figure 1: Comparison of HERMES [Q^,z] proton data on tt'*' (left) and •n~ multiplicities (right) with 
the best NLO fit curves. The error bands (the area between the dot curves) correspond to uncertainty 
estimates at 68% C.L. The errors of the data are total, statistical and systematic taken in quadrature. 


of the proton data (the mean value of per point is equal to 0.83 for tt^ and 0.65 for vr” 
multiplicities) is better than that of the deuteron data (where the mean value of per 
point is equal to 0.98 for tt’*' as well as for vr” multiplicities). The quality of the £t to the 
data is illustrated in Fig. [1] (for the proton target) and Fig. |2] (for the deuteron target). 
The error bands (the area between the dot curves) correspond to uncertainty estimates 
at 68% C.L. Note that the vertical scale is linear, not logarithmic. 


TABLE I. The parameters of the NLO input FFs at = 1 GeV^ obtained 
from the best fit to the data. The parameters marked by (*) are fixed. 


Flavor 

N 

a 

/S 

7 

h 

u 

0.278 ± 0.016 

0.276 ± 0.192 

0.188 ± 0.187 

7.64 ± 1.62 

3.09 ± 0.40 

u 

0.153 ± 0.016 

0.282 ± 0.251 

1* 

9.18 ± 4.12 

3.85 ± 0.45 

g 

0.113 ± 0.005 

12.70 ± 5.64 

14.39 ± 6.35 

0* 

- 


The values for the parameters of the input FFs (I5l) obtained from the best fit to the 
data are presented in Table 1. It turned out during the fit that there was a slight preference 
for the parameter (5^ to go to the somewhat unphysical limit zero, but the value of 
as well as the values of (z) for the measured range of z, z G [0.2, 0.8], practically 
do not change for fixed values of f3u in the range [0, 2], That is why it was fixed at the 
reasonable value = 1. Also, because of the small range of the HERMES data, a 
simpler parametrization for the gluon FF Dg^{z) was used with only three parameters 
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Figure 2: Comparison of HERMES [Q^, z] deuteron data on 7r+ (left) and tt~ multiplicities (right) with 
the best NLO fit curves. The error bands (the area between the dot curves) correspond to uncertainty 
estimates at 68% C.L. The errors of the data are total, statistical and systematic taken in quadrature. 


and 7 g = 0 . 

The extracted pion FFs from the £t to HERMES [Q‘^, z] data on pion multiplicities 
are presented in Fig. 3 along with their error bands corresponding to the uncertainty 
estimates at 68 % C.L, and compared to those determined recently by DSEHS from their 
global analysis [9] which also made use of the HERMES z] data. In Fig. 3 the 
error band for the gluon fragmentation function corresponding to = 1 (the black 
shaded band) is also presented. The corresponding error bands for the other FFs are not 
presented because they are very narrow and practically not visible. The fragmentation 
functions are plotted for the mean value of for the HERMES data, = 2.5 GeV^, 
and for the measured z region [0.2-0. 8 ]. One can see from Fig. [3] that our (LSS) pion FFs 
D'!^^ [z) and [z) are close to those of DSEHS (solid curves). In the DSEHS analysis the 
equality (|3]) for the fragmentation of the strange and u quarks into pion is not assumed. 
As a result, the extracted fragmentation functions for the strange quark, DJ^(z), differ 
a little in the z range 0.2 < 2 ; < 0.35. The main difference between the extracted FFs is 
for the gluons. This is not unexpected bearing in mind that an accurate determination 
of the qluon fragmentation function requires data covering a large range in and that 
for the semi-inclusive DIS processes the range for the HERMES data is small: 

1.1 < Q 2 < 7.4 GeV^. 

We have tried to get a feeling for the dependence of the results on the unpolarized 
PDFs used in the analysis, and End that when the MRST’02 set is used instead of the 
NLO MSTW’08 the description of the data is slightly worse, with a value of x^/d.o.f equal 
to 1.00 (0.94 for MSTW’08 PDFs). In Fig. IHwe illustrate the sensitivity of the extracted 
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Figure 3: Comparison between our pion FFs at = 2.5 GeV^ along with the uncertainty estimates at 
68% C.L. (the area between the dashed curves) and those of DSEHS (solid curves). For the gluon frag¬ 
mentation function the uncertainty corresponding to = 1 (the black shaded band) is also presented. 


pion FFs to the use of different sets of NLO unpolarized PDFs, in our case MWST’08 
and MRST’02. The corresponding FFs (z) and (z) are not shown because the 
differences between them are so small that they are not visible. Instead, for them, the 
error bands corresponding to the uncertainty estimates at 68% C.L. for D)[^(MRST'08) 
(the black solid curves) and (z) (the dashed curves), respectively, are plotted in Fig. 
4(left), and compared with the differences and , short dash and dash dot dot 

curves, respectively, where 

ad;; = D;t(MRST'02) - D;t(MSTW'08). (6) 

Also for the gluons the difference is very small, but is at least visible, as shown in Fig. 
4(right). Note that because of the large uncertainty in the determination of the gluon 
FF in Fig. 4(right) only the error band corresponding to = 1 is presented. As seen 
from Fig. 01 the central values of the fragmentation functions corresponding to the use 
of MRST’02 PDFs he entirely within the error bands for FFs corresponding to the use 
of MSTW’08 set of PDFs. The fact that a choice of PDFs other than the MSTW’08 set 
does not substantially alter the results of the global £t was mentioned also in [9]. Thus, 
to summarize, the extraction of the FFs is weakly dependent on the choice of unpolarized 
PDFs. 
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Figure 4: Sensitivity of the extracted favored and unfavored FFs (left), and gluon FF (right) to the 
choice of the set of unpolarized PDFs (see the text). Note the extremely small scale of the vertical axis 
in Fig. 4(left). 


Using the extracted FFs from the HERMES data on multiplicities in the [Q‘^,z] pre¬ 
sentation we have calculated the multiplicities at the kinematic points for the data in the 
[x, z] presentation. The obtained value for is huge, 2187.8 for 144 experimental points 
(recall that the corresponding value of for the [Q'^,z] data is 123.95). The results 
are shown in Fig. [5] for the proton and in Fig. [6] for the deuteron target. The theoret¬ 
ical multiplicities are presented along with their uncertainty estimates corresponding to 
68% C.L. As seen from the hgures, the discrepancy is very large for both the proton and 
deuteron targets for the hrst two z-bins [0.2-0.3] and [0.3-0.4], as well as at lowest x, for all 
z-bins. In our opinion such a signihcant discrepancy is totally unphysical. In an attempt 
to understand this we have tried to £t the HERMES [x, z] data directly and found that 
we cannot obtain a £t with a reasonable using different input parametrizations for the 
fragmentation functions. In addition, for some of the parameters we obtain values in the 
non-physical region. It is clear that the trend of the data in the small x region is different 
not only from that of the QCD predictions in this region, but also from the rest of the 
data points in each z bin. Consequently we decided to perform a NLO QCD £t to [x, z] 
data after removing the three lowest x data points for every z bin. The total number of 
removed data points for and tt~ multiplicities is 48 for which the contribution to 
above is 1470 (30.6 per point). In the £t to the rest of the data (96 points; we will refer to 
this data set as the ’’cut” [x, z] data) we have used for the input FFs the parametrization 
given in Eq. ([5]). Not unexpectedly it turned out that the input parameters for the gluon 
FF can not be hxed well from the £t, so for them we have used the parameters obtained 
from the £t to the [Q‘^,z] data (see Table I). The following value for x^/d.o.f, y^/d.o.f 
= 179.37/87 = 2.06 for 96 experimental points and 9 free parameters, is achieved in the 
£t. The results of the best fit are shown in Fig. 5 for a proton target and Fig. 6 for the 
deuteron one (solid curves). Their continuation to the low x region where the data points 
were removed from the fit, is indicated by the dashed curves. The quality of this fit is 
illustrated in Table H, and compared to the quality achieved in the fit to the [Q^, z] data. 
It follows from the values, presented in Table H, that the description of the [Q^, z] data 
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Figure 5: Comparison of HERMES [x,z\ proton data on 7r+ (left) and 7r“ multiplicities (right) with 
the multiplicities at the same kinematic points calculated by our EEs extracted from HERMES [Q'^^z] 
data (dot curves with the bands corresponding to uncertainty estimates at 68% C.L). The errors of the 
data are total, statistical and systematic taken in quadrature. The solid curves correspond to the best fit 
to the cut [x, zj data (see the text). 


is much better than that of the [x, z] data even after removing a third of the data points. 

The extracted pion favored and unfavored FFs from the fit to the HERMES cut [x, z] 
data on pion multiplicities are presented in Fig. 7 and compared to those determined 
from the fit to data, for which the error bands corresponding to the uncertainty 

estimates at 68% C.L. are also presented. Recall that the gluon FF is the same for both 
the representations of the data and it is shown in Fig. 3. As seen from Fig. 7, the 
central values of the favored pion FF (solid curve) extracted from the cut [x, z] data are 
systematically smaller than those extracted from [Q‘^,z] data, and in the z region [0.2, 
0.4] the corresponding curve lies outside the error band. The central values of unfavored 
pion FF extracted from the cut [x, z] data are also systematically smaller then those 
extracted from data, however, the corresponding curve lies within the error band. 

It is important to mention, however, that from the calculation of the multiplicities at the 
kinematic points for the data in the [Q^, z] presentation, using the extracted FFs from the 
fit to the cut [x, z] data, we obtain for the value 665.2 which is more than five times 
larger than the value 123.95 achieved in the direct ht to the [Q‘^,z] data. 
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Figure 6: As in Fig. 5 but for a deuteron target. 


TABLE II. per point values for the pion multiplicities 
obtained from the hts to [Q^, z] and the cut [x, z] data. 



fit 

[x, z] fit 

P 

0.83 

1.94 

Mf 

0.65 

1.58 

Mf 

0.98 

1.63 

Mf 

0.98 

2.33 


Note that in all our NLO QCD calculations of the [x, z] pion multiplicities we have 
used for x and their mean values < x > and < > as given in the HERMES 

data tables |5]. We have checked that in NLO QCD the pion multiplicities calculated 
at the average kinematics {< x >,< >} coincide extremely closely (to better than 

1%) with the average multiplicities calculated using the expression, Eq. (1), given in the 
recent HERMES paper ra as applied to the NLO semi-inclusive and DIS cross sections 
(see the remark m)- This fact is very important because it means that the huge time 
consuming the computer calculations involved in using the above mentioned expression 
in fitting the data on the average multiplicities can be significantly reduced if the NLO 
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QCD multiplicities are calculated at the corresponding mean values < x > and < >. 

Finally we would like to underline that our NLO QCD analysis of the HERMES 
[x,z\ data supports the assertion of Stolarski |T5], based on a LO QCD analysis, that 
the increase of magnitude of the HERMES pion multiplicitity sum as x decreases in the 
region x < 0.1, is difficult to reconcile with perturbative QCD. While in [TJ] the argument 
is presented for a deuteron target, our observation is that it holds for both proton and 
deuteron targets. 



Figure 7: Comparison between favored (left) and unfavored (right) fragmentation func¬ 
tions extracted from [Q^, z] (dot curves) and the cut [x, z] HERMES data on multiplicities 
(solid curves). The dashed curves mark the error bands for FFs determined from the 
[Q^, z] data. 


4 Summary 

The publication by HERMES of the final version of their data on pion multiplicities 
on protons and deuterons has profound implications for our understanding of the pion 
fragmentation functions. 

1) The fact that the final data are significantly different from the preliminary data 
means that the oft utilized DSS FFs |2], which were based on the preliminary data, are 
incorrect. 

2) We have studied the two-dimensional projections of the hnal HERMES data, the 
so-called [x,z\ and [Q‘^,z] formats, presented by the HERMES group. 

a) With the pion FFs, parametrized in a standard way, and respecting isospin invariance, 
we have found an excellent fit to the [Q^, z] presentation of the data and extracted a new 
set of NLO pion FFs. Except for the gluon fragmentation function, our new pion FFs 
are very similar indeed to those obtained recently by the DSEHS group [9] using, in their 
global analysis, the [Q‘^,z] HERMES data. 

b) On the contrary, no reasonable NLO QCD fit could be achieved to the [x, z] 
presentation of the data. We have found that an adequate ht to the [x, z] data is only 
possible if we cut points with x < 0.075 from the data which means that a third of the data 
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points is removed. However, even with these cuts, the quality of the description of the 
[Q^, z] data is much better than that achieved for the cut [x, z] data. While the extracted 
unfavored pion FF lies within the error band corresponding to the unfavored pion FF 
extracted from the £t to the [Q^,z] data, the favored pion FF[a;, 2 :] is systematically 
smaller than favored FF[(5^, z] and is outside of its error band in the region 0.2 < < 0.4. 

3) We have found that the trend of the data in the HERMES [x, z] presentation of 
their data, where the magnitude of the pion multiplicitities in the region a; < 0.1 increases 
as X decreases, is totally at variance with the trend of the NLO QCD predictions. This 
suggests that possibly there is a problem with the HERMES [x, z] presentation of their 
data, and emphasizes the need for new data on the hadron multiplicities. We thus await 
with great interest the publication of the final COMPASS data on the pion multiplicities. 
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